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We consider a cold gas of non-interacting fermions in a two dimensional harmonic trap with two 
different trapping frequencies uj^ ^ Uy, and discuss the effect of rotation on the density profile. De- 
pending on the rotation frequency and the trap anisotropy Uy/uix, the density profile assumes two 
qualitatively different shapes. For small anisotropy {ujy/ujx «C ^f + A^'^ /u)%), the density consists 
of elliptical plateaus of constant density, corresponding to Landau levels and is well described by a 
two dimensional local density approximation. For large anisotropy {ujy/uJx -y/l + 4S72/lj|), the 
4_> density profile is Gaussian in the strong confining direction and semicircular with prominent Friedel 

O oscillations in the weak direction. In this regime, a one dimensional local density approximation is 

well suited to describe the system. The crossover between the two regimes is smooth where the step 
structure between the Landau level edges turn into Friedel oscillations. Increasing the temperature 
causes the step structure or the Friedel oscillations to wash out leaving a Boltzmann gas density 
profile. 
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I. INTRODUCTION 
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^ Rotating ultracold gases have received a lot of attention after the initial experimental demonstration of vortices 

and vortex lattices in a rotating Bose Einstein condensate [TJ [3]. Beyond vortex physics these experiments hold 
I ; promise to create clean, controllable environments to investigate interesting quantum phases brought upon by the 
unusually high effective magnetic fields created by rotation. There are both theoretical proposals and experimental 
^ efforts to investigate vortex lattice melting [U [5], Fractional Quantum Hall states [BJ|71IHl[n] and BEC-BCS crossover 
I for rotating gases (TUJ EH H US] . 

^ The usual procedure for rotating gas experiments start by pumping large amount of angular momentum into the 

^ cloud by various rotation methods [21 13] . Later, confining the gas into trap that is as isotropic as possible to preserve 
Q angular momentum, the system settles into a ground state that is rotating at some equilibrium frequency. Another 
I— I alternative route, which is used in at least one experiment f 141, is to create a rotating anisotropic trap with some fixed 
rotation frequency and wait for the system to reach the ground state at this fixed rotation frequency. This procedure 
^ allows for one more control parameter, the trap anisotropy, and is theoretically interesting as it allows the system to 
I smoothly go from a two dimensional regime for small anisotropy to a one dimensional regime for large anisotropy. 
\^ The single particle Hamiltonian for a rotating anisotropic trap is exactly diagonalizable, where the eigenenergies can 

be found by a Bogoluibov transformation [T5l [TB]. However, a closed form expression for the real space wavefunctions 
•/^ has not been given until recently 17]. The investigations of rotating anisotropic trap systems have focused either on 
the properties of vortices, or dynamics in the extremely anisotropic limit |18l [TOl ED] . All of these studies concern the 
y—^ behavior of interacting Bosons in an anisotropic rotating trap. The interaction strength creates another energy scale 
which complicates the crossover between the one dimensional and two dimensional regimes of the anisotropic trap. 
In this study, we consider the density profile for a gas of non-interacting fermions in a rotating anisotropic trap. 
^ Because of the Pauli suppression of s- wave scattering, spin polarized ultracold fermi gases are naturally non-interacting. 
The absence of another energy scale related to interactions simplifies the physics considerably and reveals how the 
anisotropic trap smoothly connects the one dimensional regime to two dimensional regime. 
^ The density profile for non-interacting fermions are easily calculated numerically by summing the densities of 

individual eigenstates. For the rotating anisotropic trap the wavefunctions consist of a exponential part multiplying a 
Hermite polynomial of a complex coordinate [17J. The recursion relation for Hermite polynomials allows for fast and 
accurate numerical calculation of the density, enabling us to calculate the exact density for thousands of particles for 
any rotation frequency or anisotropy. 
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The numerically calculated density profiles clearly show two distinct regimes. For small anisotropy the density 
profile consists of a ziggurat like structure of elliptical density plateaus. The density inside a plateau is an integer 
times ^ where 7 is a function of trap frequencies and rotation. Each plateau is clearly identified to belong to 
a Landau Level, and this regime can be understood as a small deformation of the quantized density profile in an 
isotropic rotating trap which was first discusses by Ho and Ciobanu[.21 . As the trap anisotropy is increased, the 
plateau structure is replaced by a profile that is Gaussian in the strongly confined direction and a sum of semicircles 
in the weakly confined direction. Friedel oscillations of the density become prominent as this switching happens, 
and for very large anisotropy the Landau Level quantization is replaced by quantization of the wavefunction along 
the transverse, i.e strongly confined, direction. The density profiles give a striking example for smooth connection 
between the two dimensional and one dimensional physics. 

The paper is organized as follows: We review the single particle energies and wave functions in the next section and 
give the expressions for the wavefunctions in the lowest five Landau levels. In section III we show that the crossover 
from two to one dimension is already evident in the energy spectra and define the two regimes quantitatively. Section 
rV focuses on the numerical density profiles and the description of the properties in both regimes. We also investigate 
the density profile analytically, and discuss one dimensional and two dimensional local density approximation methods 
along with Friedel oscillations. Section V contains the numerical results at finite temperature. We conclude in section 
VI by briefly discussing the experimental consequences of our results and indicating directions for future investigations. 



II. SINGLE PARTICLE PROBLEM 



We start by considering the single particle physics in a two dimensional rotating anisotropic trap. In this section, 
we review the results obtained by Fetter, who obtained expressions for the wavefunctions in the Lowest Landau level 
in reference |17j . 

Our discussion throughout the paper will be limited to behavior in two dimensions. With this limitation, we are able 
to focus on the physics of switching between one and two dimensional regimes without the complications introduced 
by the motion in the third dimension. Our results can be extended to three dimensional systems quite easily if there 
is strong confinement in the third dimension or if the potential in the third dimension is slowly varying [3T]. Both of 
these limits are routinely realized in ultracold atom experiments. 

The Hamiltonian for a single particle in a two dimensional anisotropic harmonic trap can be written in the rotating 
frame as 

^ - ^ (P' + pI) + \mujIx^ + \Mu:y - ni,. (1) 

Here M is the mass of the particle, the angular momentum in the z direction Lz is given as — xpy — ypx, and U, is 
the rotation frequency. The trapping frequencies along the x and y directions are w^, Wy, respectively. Without loss of 
generality we take ujy > oJx, and refer to the y and x directions as the strongly confined and weakly confined directions 
throughout the paper. The stability of the trapped system depends on the rotation frequency, which requires the the 
rotation frequency to be smaller than the weak confining frequency < w^. 

It is clear from the single particle Hamiltonian that the system is essentially described by two dimensionless pa- 
rameters, the anisotropy 



and the scaled rotation frequency 



(i, = > 1, (2) 



n = — <i. (3) 



All the relevant quantit ies can be non-dimensionalized by scaling all the lengths by the oscillator length in the 
X direction = \/ h/mujx, frequencies by ujx and energies by hujx- In the recent literature on rotating anisotropic 
traps |T6j [20] different sets of non-dimensional parameters have been used to describe the system, unnecessarily com- 
plicating the interpretation of the results. We avoid increasing this list of dimensionless quantities by following the 
notation of Fetter|17j closely. Our numerical calculations are done in dimensionless units, and we will explicitly state 
our scaling where appropriate. 

The eigenenergies for the Hamiltonian Eq.Q can be obtained by a direct diagonalization, or two successive Bogolui- 
bov transformations [15, tl6j. The eigenstates are labelled by two integers n > and m > 0, and the corresponding 
energy is 



In r ^ ' 



(4) 
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where 



W4 



(5) 



J). A schematic plot of the energy levels is given in figure 1, where one can observe that if the 
rotation is rapid U, ^ ojx, then u>+ <C OJ- and the levels with the same index n are almost degenerate. Using the 
analogy between a rotating system and a system under a magnetic field, we identify the levels with the same n index 
to belong to a Landau level, and call n the Landau level index. 

While the system is controlled by the two dimensionless parameters in Eqs. ([2]), (|3]), a number of other derived 
parameters give insight into the dynamics and allow one to express results more concisely. The anisotropy of single 
particle orbits are controlled by two dimensionless parameters 
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both of which are between and 1. We show in the following sections that a particular combination of these parameters 
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controls the switching between the one and two dimensional regimes. The density of particles in a filled Landau level 
(see Section IV) is related to the frequency 
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Finally we define the two 'stretched' complex coordinates which control motion in Landau levels 

/2A/7/3+ x + if3^y 
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The ground state of the system has the wavefunction 
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where and ay are the widths of the Gaussian envelope 

2 i + P+P- h 



1 + l3+p. 



1 + f3+(3_ 



f3+ 



(10) 
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The wavefunctions in the lowest Landau level (LLL) are found by applying the relevant raising operator to be 



(11) 



(t>mo{.x,y) = 



1 



m/2 



(1^00 ix,y) 



(12) 



where Hm is the to*'' Hermite polynomial. 

Another raising operator allows one to generate the wavefunctions in the higher Landau Levels, for example the 
states in the first excited Landau level have the wavefunction 



0ml (a;, y) 
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with p = y/0lp^/{l + l3+f3_). We give the wavefunctions for the six lowest Landau levels in the appendix, and use 
them in the numerical calculations in section III. 
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FIG. 1: Energy vs. number of states in each Landau levels. / is the highest number of states within LL's before the next LL is 
excited. 



The wavefunctions given above are remarkable, as they continuously connect the usual ID harmonic oscillator 
wavefunctions which have the Hermite polynomial form to 2D Landau level wavefunctions which are (for the lowest 
Landau level) analytic functions of the complex coordinate z = x + iy. In the isotropic limit the m*'' wavefunction 
has a m*'' order zero at the origin. As the anisotropy is turned on the m*'' order zero immediately breaks up into m 
first order zeros all of which are on the x axis, as the roots of Hermite polynomials are always real|16j. Increasing 
anisotropy separates the roots from each other and forces the Gaussian envelope to be more and more anisotropic. In 
the next section we analyze the one and two dimensional regimes by concentrating on a single quantity that is found 
from the eigenvalues: the number of energy states in the LLL that have less energy than the first excited Landau 
Level. 



III. ONE AND TWO DIMENSIONAL REGIMES 



The switching between the two dimensional and one dimensional regimes is already apparent in the energy spectrum 
given m Eq.Q. It is institutive to consider two limiting cases. The energies of eigenstates of a rotating isotropic trap 
given as 
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1 



The energies for a non-rotating anisotropic trap are 



Er, 



+ -)UJ^ + {Uy + -}ujy 



(14) 



(15) 



Both of these expressions have the same structure, as they are labeled by two non-negative integers, and this structure 
can be graphically represented as in figure 1. The rotating anisotropic trap system interpolates between these two 
limits, and a good identifier for this interpolation is the number of lowest Landau Level states that have lower energy 
than any state in the higher Landau levels. We define 



I : 



UJ+'' 



(16) 



which is graphically shown in figure 1. 

The behavior of / as a function of LOy/uJx at fixed rotation frequency ^/lUx clearly shows the two regimes (See 
Fig. (2)). At zero anisotropy, I has its isotropic value (1 +Ct)/{1 — Cl). When the anisotropy is introduced, this value 
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FIG. 2: Solid line shows the behavior of 1, number of the stats within two LL's as a function of the anisotropy of the system 
uiy/Q,x- Dashed lines are the asymptotic results for one and two dimensional limits. Hero ujyc = (1 + AQ,/ijjxY'^ is the critical 
anisotropy where 1 is minimum. The minimum number of states is given by Ic = ( i-n/wa, ■^^'^^ 



quickly decreases as {ijy — 1) and obtains its minimum value of Zmin = ^^8/(1 — Q) at = ^/Y+~4D?. As the 

anisotropy is increased further, I increases Hnearly with ujy with slope 1 — This Hnear increase is exactly 
what is expected in the one dimensional regime as here / would be the ratio of the effective frequencies in the y and 
X directions, 



/ W2 -0,2 - 

The degeneracy of the Landau levels in the two regimes originate from two different physical effects. When the 

anisotropy is small, and the system is in the two dimensional regime, Coriolis force induced by the rotation partially 
freezes the kinetic energy of the particles and causes the high degeneracy within the Landau levels. For high anisotropy 
one can think of the energy spectrum in terms of subband quantization in a narrow channel, where each 'Landau 
Level' actually is formed by stat(^s that have the same wavefunction in the tightly confined direction. 

The form of the wavefunctions also support this identification based on I. For low anisotropy the wavefunctions 
peak around an ellipse in the x — y plane, while in the one dimensional regime the density in the x direction shows 
prominent oscillations, similar to a one dimensional harmonic oscillator wavefunctions. In the next section, we show 
that the density of a Fermi gas trapped in a rotating anisotropic trap has markedly different behavior in these two 
regimes. 



IV. DENSITY PROFILES 



We now consider the density of a gas of A'' non-interacting spinless fermions in a rotating anisotropic trap. As 

the real space wavefunctions are known the calculation of density reduces to a sum over the absolute squares of the 
wavefunctions of all filled eigenstates. Hence the density is given by 

Pi^^y) =^\4'nm{x,y)f 0{lJ - Enm), (18) 
nm 

where fj, is the chemical potential and the energies Enm and corresponding wavefunctions ^nm are given in the previous 
sections. The chemical potential determined by the number of particles N in the system with the constraint 



dxdyp{x, y) = N. 



(19) 
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The LLL wavefunctions given in sec II are composed of an exponential factor multiplying a Hermite polynomial, and 
wavunctions in the first few Landau levels given in the appendix also have a similar structure. Hermite polynomials, 
like other orthogonal polynomial satisfy a recursion relation [22] 



Hn+i{x) = 2xHn{x) - 2nHn-i{x), 



(20) 



which enables rapid and accurate numerical evaluation of the wavefunctions. Because of the recursive structure, the 
density sum Eq.( 18 ) can also be calculated numerically at very low computational cost. We calculated density profiles 
for up to = 2000 particles, for parameter values such that all the particles reside in the lowest six Landau Levels. 

In figure 3, we display the density profiles for N = 1000 particles with fixed rotation frequency ^/lOx — 0.999, for 
different values of anisotropy. These parameters demonstrate the general behavior of the density profile when the 
particles are distributed over the first few Landau levels. For small anisotropy the density profile mainly consists of 
steps, with small regions of switching between them. Each plateau is clearly linked with a Landau level, and defines 
an elliptical area of almost constant density. As shown below, the density contribution from a filled Landau level is 
p = Thus, the steps have density equal to an integer times For the parameters in Fig. (3a) only the lowest 
three Landau levels have particles and corresponding steps can be identified clearly. The switching between the n'^ 
step and (n — l)"* steps happens with ti — 1 oscillations in density, which can also be seen in the figure. As the 
anisotropy is increased the density oscillations between the Landau level steps become more prominent especially in 
the weak trapping direction. When the anisotropy is within the critical region of switching between two and one 
dimensional behavior, the steps and the Landau level structure become smeared out. Instead, one can observe the 
density oscillations becoming prominent throughout the cloud. These oscillations are expected for Fermions, as a 
sharp Fermi surface cutoff in momentum space results in oscillations in real space, known as Friedel oscillations |23j. 
For larger values of anisotropy, the density profile assumes a Gaussian shape in the strong confining direction and the 
familiar semicircular shape of one dimensional trapped fermions in the weak confining direction. Each Landau Level, 
or rather sub-band in this regime, adds a new semicircle in this direction, thus it is quite easy to observe how many 
sub-bands are occupied form the density profile (See Fig.(3f)). 

In figure 4 we display the density profiles with N=200, a low enough particle number so that all the particles remain 
in the LLL for any value of anisotropy. In this figure, it is easier to discern the character of the change between the 
two and one dimensional regimes, as one can observe how the flat plateau is replaced by a semicircular shape with 
Friedel oscillations. Again, the switching happens around the critical value of anisotropy determined from I. 

To gain a better understanding of the density profile, we investigate the density sum Eq.(18l in more detail. Let's 



concentrate on the case for which all the atoms reside in the LLL, for N particles we can write the sum as 



N 



m=0 

^ to! V2/ V2c v2c TTa^ay 



(21) 



m— V V y 

In the limit N oc this sum can be evaluated using the generating function for Hermite polynomials. We define 



(22) 



We use the generating function for Hermite polynomials 



-t^+2tz 
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Hyn{z) 



(23) 



but regard t = |t|e*^ as a complex variable. Then the infinite sum can be written as 

r2^ ,oo 



S{z,z) = 
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The resulting Gaussian integral can be evaluated to yield 



5(z, z) 
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FIG. 3: (color online) Density profiles for fermions at different anisotropy values of the system and fixed number of the particles 
N=1000, and fixed rotation frequency Q/uJx = 0.999. The Friedel oscillations are observed in density profiles of anisotropic 
cases. The number of LL's that fermions fill at each anisotropy can also be determined by counting the number of plateaus in 
the density profiles. 



which gives the constant density 



p{x,y) = Po = 



(26) 



We see that 7 is the parameter that controls the density of states per Landau level in an anisotropic trap. In the 
numerical calculations we see that the steps in the density profile appear at integer multiples of this value. In the 
two dimensional regime, when the anisotropy is small, the density at the center of the trap quickly reaches this value. 
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FIG. 4; (color online) Density profiles of fermions at different anisotropy values of the system and fixed number of the particles 
N=200, and fixed rotation parameter Q,/uxa = 0.9999, when all the fermions are settled at the LLL. The Friedel oscillations 
are also observed for density profiles of anisotropic cases. 

even for small number of particles. This is because the contribution of higher wavefunctions to the density near the 
center is small for low anisotropy. 

We must remark that the above result is correct for any anisotropy, i.e. if one puts infinitely many particles into the 
LLL one would get a constant density. However as the particle number is increased the cloud goes through two very 
different paths to this constant density profile. In the two dimensional regime, the sum of the density of N particles 
defines an elliptical region of constant density which grows in size as N is increased. In the one dimensional regime, 
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density is very anisotropic and shows prominent Friedel oscillations. The frequency of these oscillations decreases as 
N is increased and the density becomes smooth very slowly with ^ 

The presence of Friedel oscillations in a trapped Fermion system is not unexpected, as the physical reason for their 
presence is the sharp cutoff at Fermi energy. When one is in the two dimensional Landau Level regime the oscillations 
are suppressed as locally the maximum allowed density within a Landau level is reached. One can investigate the 
Friedel oscillations by a simple scaling argument as follows. First let us define the finite sum over Hcrmite polynomials 
as 

SnMc)^}_^ — c, (27) 

with z = a-\-ih. Now because of the infinite sum result we have obtained above we expect two such sums to be related 
as 



S„(»,..).j'i^S„(Jll±£lg„Jii^MO, (28) 



if A'' ^ 1 is sufficiently large. Now ci can be chosen arbitrarily close to one, 

N 



SN{a,b,c,) = ^ — c. 



H„{z)H,,iz) 



n=0 
N 

^ ^ 2"n! 

ri=0 

= SN{a,b,l), (29) 
provided that (1 — ci)N ^ 1. This last sum can be evaluated exactly using Christoffel-Darboux formula 

ri=0 

^ Hn+i{z)Hm{z) - Hn+i{z)Hn{z) 

2^+iiV!(z-z) ■ ^ ' 

This exact evaluation and the scaling assumption readily yields a formula for the wavevector of the Friedel oscillations 
near the center of the trap 

p{x,y^ 0) = po (^1 + i-lf ^cosikpx)^ , (31) 

with 

kpl^ = ^/m/LJ^^/L0++L0^P+l3~. (32) 

As expected, this scaling result overestimates the amplitude of the Friedel oscillations, especially in the two dimensional 
regime. The suppression of oscillations due to limited local density of states is not captured by the scaling assumption. 
However, scaling correctly describes the wavelength of the oscillations well into the regime where they are too small 
to observe within our numerical precision. 

When the external potential is smooth, a very useful tool to describe the density profile is local density approxi- 
mation. For the case of a rotating anisotropic lattice we need two different versions of LDA to describe the different 
behavior in the two regimes. 

For the two dimensional regime the effective potential in both strong and weak confining directions are smooth, 
when compared with the magnetic length imposed by rotation. Thus, it becomes permissible to treat the system 
locally as a homogenous system with the Landau Level structure. The local density at any point of the gas will be an 
integer multiple of the density of a filled Landau level mfi / h. Thus the density profile within the 2D LDA consists of 
a sequence of elliptical plateaus, which matches well with the numerically calculated density profiles in the 2D regime. 
Thus 2D LDA is successful in describing the density profile except for the switching regions between the Landau level 
steps. A typical plot comparing the result of LDA with the calculated density profile is given in Fig. (6). 
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FIG. 5: (a) Density for N=200 fcrmions at n/t^^ = 0.9999, and = 2.2. (b) Density for N=300 ferrnions at fi/wx = 0.999, 

and ujy/Q.x = 2.2. The insets in both figures highlight the Friedel osciUation in density profiles of the ferrnions. 




FIG. 6: (a) Density of N=200 ferrnions obtained by direct numerical calculations (solid line), and one dimensional LDA (dashed 
line) at fi/wa, = 0.9999, and ijjyj0.x = 3.5. (b) Density of N=325 fermions obtained by direct numerical calculations (solid line), 
and two dimensional LDA (dashed line) at 0./ijJx = 0.99, and u>y/Q,x = 1-2. 



As the anisotropy is increased and the gas enters the crossover regime the step structure in the density profile is no 
longer observed. This regime is not well described by a two dimensional LDA, as the potential in the strong confining 
direction is no longer smooth on the scale of magnetic length. However if the anisotropy is increased further, the 
oscillator length in the strong confining direction becomes much smaller than the magnetic length and another LDA 
approach becomes feasible. This is exactly the one dimensional regime described in Sec. IIL 

In the one dimensional regime, the density profile in the strong confining direction is almost the same for all 
wavefunctions in the same Landau level. Thus we can treat the system within a one dimensional LDA and obtain the 
usual semicircular profile for the one dimensional fermions 



where the radius of the cloud is found as 



p{x,y = Q)=p{Q,Q)\jl-j^, (33) 
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FIG. 7: (a) Density profile for ferinioiis at different temperatures and at Q.IuJx = 0.99, and uiy/Vlx = 3.5. (b) Density profile 
for fermions at different temperatures and at Q,/uix = 0.99, and ujyj^x = 1.01. The oscillations (a) and the layered aspect of 
the profiles (b) are eliminated by temperature. Here /3 = is the normalized inverse temperature. 

As can be observed in Fig. (6), this approximation describes the density profile in the one dimensional regime rea- 
sonably. The agreement increases with increasing anisotropy, but Friedel oscillations can not be captured within the 
local density approximation. The agreement of one dimensional LDA with the numerically calculated profiles further 
support our interpretation of the very anisotropic Landau levels as sub-bands formed by quantization in the strong 
confinement direction. 



V. TEMPERATURE EFFECTS 

The density profiles given in the previous sections were all calculated at zero temperature. At finite temperature 
the density profile can be calculated by including the fermi distribution function as 

P{x,y) = ,,p^(^^^_^) + i . (35) 

where /3 = is the inverse temperature. 

We calculated density profiles for finite temperatures using the same numerical procedure with the zero temperature 

case. The temperatures and chemical potentials were chosen to ensure that the occupation of the states higher than 
the sixth Landau level is less than ea;p[— 5] ~ 0.0067. For the two dimensional regime the step structure of the density 
profile is expected to be smeared out when the temperature is enough to close the gap between the Landau levels 
ksT ^ huj^. The calculated density profiles show that the step structure becomes indiscernable below this value, 
as the transition regions between the Landau level steps dominate the profile. For larger temperatures the profile 
assumes the gaussian profile as would be expected from a Boltzmann gas. 

For the one dimensional regime, the first effect of the finite temperature is to smear out the Friedel oscillations. The 
relevant energy scale for Friedel oscillations is the separation of energy levels within a Landau level, and as soon 
as this temperature is reached the density profile becomes smooth, while still obeying the semicircular LDA result. 
For higher temperatures, once again the density profile becomes Gaussian as would be expected from a Boltzmann 
gas. The effects of temperature in both regimes is plotted for typical density profiles in figure 7. 

VI. CONCLUSION 

We show that the density profile of a non-interacting Fermi gas in a rotating anisotropic trap has two distinct regimes 
depending on the rotation rate and anisotropy. Through numerical calculation and a local density approximation the 

distinct behavior in the one and two dimensional regimes is observed. For small anisotropy (cJj^/cJj. ^ -^/l + 4i72/ti;^), 
the density consists of elliptical plateaus of constant density, corresponding to Landau levels and is well described by 
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a two dimensional local density approximation. For large anisotropy [ujy/LOx ^ a/1 + /lo'^)^ the density profile is 
Gaussian in the strong confining direction and semicircular with prominent Friedel oscillations in the weak direction. In 
this regime, a one dimensional local density approximation is well suited to describe the system. The crossover between 
the two regimes is smooth where the step structure between the Landau level edges turn into Friedel oscillations. 
Increasing temperature smears out the step structure in the two dimensional regime, and similarly smoothes out the 
Friedel oscillation in the one dimensional regime. 

The exact wavefunctions given for the rotating anisotropic trap provide a smooth one to one mapping between 
the states in a 2D Landau level and the states in a one dimensional harmonic trap. We expect this mapping to be 
useful for study of interacting Bose and Fermi systems, going beyond the non-interacting gas in this paper. It would 
be interesting to see what kind of states are obtained for an extremely elongated trap if this mapping is used on 
correlated states such as the fractional quantum Hall states in the two dimensional regime. Similarly, certain strongly 
interacting models in one dimension may accept continuation into two dimension, extending the validity of some exact 
solutions. 

Experiments on rotating ultracold gases have so far focused on isotropic traps, as they mostly rely on the isotropic 
trap to conserve angular momentum given to the system before equilibrium is established. Nevertheless a rotating 
anisotropic trap has been demonstrated [1 J. A wide array of experiments have displayed the feasibility of designing 
traps that are strongly confining in one or two dimensions to probe low dimensional physics of quantum gases. All 
these successes with trap design lead us to conclude that a the trap potential discussed in this paper is within the 
current experimental capabilities. 

A major difficulty in Fermion experiments is creating low enough temperatures to probe quantum mechanical nature 
of the gas|25j |26l [27l [28j [29j [30]. However as the system considered in this paper is non-interacting, the temperatures 
needed to observe qualitative features such as density steps is well above the currently obtained lowest temperatures. 

The main method employed in measuring density profiles is expansion imaging of the gas, and this method can 
be used to determine the main features of the density profile such as Landau level step structure or the semicircular 
density profile in the one dimensional regime. More subtle features, such as Friedel oscillations would require more 
complicated probes of density such as Bragg spectroscopy [ST . 

Finally, we remark that the Fermi gas in a rotating anisotropic trap is system that may elucidate the connection 
between one dimensional physics with dynamics within a Landau level. It would be interesting to extend the study 
in this paper to interacting fermion systems. 



VII. APPENDIX 



The ground state wave function of the system is written as 
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(36) 



and the excited states are obtained by applying the relevant raising operators defined by Fetter[T7]on ground state 
wavefunction. Operator a+gives the The excited sates wave functions in each LL, and another one i.e. a_ gives the 
higher LL states wavefunctions. The wave functions in LLL are obtained to be 



where 



^/2c 



The wavefunctions in higher excited Landau levels also are found to be 

1 



Vm{x,y) 



^_P„.(e+) -2mipP,„_i(e+) 



(37) 



(38) 



(39) 



'Pni2{x,y) = -^=(poo{x,y)l^(^£,'i - cjPrniM 

- 4m(TO-l)p2p„_2(e+)l, 



(40) 



- 6imp(e^ -c)p„_i(e+) 

- 12m(m-l)p2^_P„_2(^+) 

+ 8im(m - l)(m - 2)p'^ Pm-3{£,+) 



- 8m^p(e--3cC_)p„_l(e+) 

- 24m(m - - c)p„_2(C+) 
+ 32im(m - l)(m - 2)p3^_P„_3(^+) 

+ 16m(m-l)(m-2)(m-3)p4p„_4(e+)| 



i=(/Poo(x,y)|(e - lOce + 15c2^_)p„(^+) 



V120TO 

- 40m(m - - 3ce_)p„_2(C+) 

+ 80im(m-l)(m-2)p3(e-c)P™-3(^+) 
+ 80m(m-l)(m-2)(m-3)p^^_P„_4(^+) 

- 32im{m - l)(m - 2)(m - 3)(m - 4)p-'^P„_5(^+) 

We suggest a general form for the wavefunctions of the excited states: 



^mn{x,y) = 
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fc=0 



^n—k 
n—k n^~^ 

{n-k)V 



(^n-fe^_ -^"(^-) ^n-fe^_^ -^™(^+) 
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